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Abstract 

One-dimensional scattering mediated by non-Hermitian Hamiltonians is studied. A 
schematic set of models is used which simulate two point interactions at a variable 
strength and distance. The feasibility of the exact construction of the amplitudes 
is achieved via the discretization of the coordinate. By direct construction it is 
shown that in all our models the probability is conserved. This feature is tentatively 
attributed to the space- and time-reflection symmetry (a.k.a. VT— symmetry) of our 
specific Hamiltonians. 
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1 Introduction 

In the absence of an external potential, the motion of a quantum particle is described 
by the kinetic-energy Hamiltonian Hq = —d 2 /dx 2 in one dimension (h = 2m = 1). 
This operator is Hermitian and, incidentally, symmetric with respect to the space and 



1 e-mail: znojil@ujf.cas.cz 



time reflection (i.e., VT — symmetric, H VT = VTH Q , cf. many relevant comments 
on such a type of symmetry in p] ) . 

In an approximation where the real line is replaced by the mere discrete lattice 
of coordinates with some sufficiently small stepsize h > 0, 



k = 0,±1, 



the role of the kinetic energy is often being played by the doubly infinite tridiagonal 
matrices 



2 -1 
-1 2 -1 
-1 2 



or H n 



-1 
-1 -1 
-1 



which differ just by a trivial shift of the energy scale. Whenever we treat V as 
the parity (Vxk = X-k) and the antilinear operator T as the time reversal (i.e., in 
our present setting, transposition plus complex conjugation), we may represent the 
product-operator symmetry of our real matrices Hq by the antidiagonal unit matrix 



(1) 



Using this definition we shall demand that also all the nontrivial, doubly infinite 
discrete Hamiltonians H = Hq + W possessing a nonvanishing interaction term W 
will be required real and VT— symmetric. 

The matrix dimension of the interaction matrix W (i.e., the "range" of the in- 
teraction) will be assumed finite. One expects that then the scattered states could 



stay asymptotically undistorted. In the mathematical terminology such an expec- 
tation means that we feel allowed to search for the solutions of the discrete and 
VT— symmetric Schrodinger equations 

(H + W)ij = Eij (2) 

complemented by the standard, undistorted boundary conditions 

I e im( ? + Re- im( ? , m < -M , 

i>m={ " (3) 

I Te imv , m>M. 

We should remind the readers that the standard re-parametrization of the energy 
E = (2 — 2 cos (f)/h 2 in terms of the real angle (p G (0, it) should be used [2J. 

Our study has been inspired by a few papers on the scattering in non-Hermitian 
scenario [3], IU [5] and, in particular, by the Jones' paper [Bj. Unfortunately, its 
author worked in the differential-equation limit h — > which made the detailed 
analysis perceivably hindered by the non-Hermiticity of the equations. In effect, the 
feasibility requirements (cf. [8]) restricted his attention to the mere VT — asymmetric 
delta-function interactions, therefore. 

In our subsequent comment [7] we facilitated the technicalities by the transition 
to the discretized eq. (J2J). Having preserved the Jones' philosophy we choose just the 
VT— asymmetric models exemplified by the "ultralocal" , two-by-two matrix example 



W (UL) 



\ a / 



such that W( UL ' VT 7^ VTW^^. Due to the discretization approximation h > 
we were able to construct the explicit formulae for the reflection and transmission 
coefficients R and T, respectively, 

R (UL ) = _^^ t (ul )= (l-a)(l-e 21 ^ A = 1 _ (1 _ a 2 )e 2i^ 



We were also able to mimic the key features of the Jones' first-order perturbation 
results by another entirely exact and compact formula 

R{UL) * + TiUL) * = l-a[l + U(a^ ^ v) _ 



l + a[l + ?7(a,^)]- 1 ' v '^ 2 (1 - a) (1 - cos 2<p) 

This formula closely resembled the very similar Jones' perturbation results [6] . Hence, 
we could also parallel his conclusion that since the predicted sum appears greater 
than 1 or less than 1 (depending on the sign of the coupling a) it cannot be given 
the usual probabilistic interpretation. One must rather assume the presence of some 
respective "unknown source" or "unknown absorber" near the origin. Thus, in the 
effective-theory manner, the mathematical non-Hermiticity of the interaction terms 
W precisely reflects the presence of certain hidden physical mechanisms which violate 
the conservation of the number of particles. 

In the context of the internal physical consistency of many non-Hermitian bound- 
state models pG| such an effective-theory physical interpretation of the scattering looks 
rather unsatisfactory. In what follows, for this reason, we shall try to re-install the 
VT— symmetry in our matrix model(s) and study the consequences. For this purpose 
we shall make use of the enhancement of the feasibility of the calculations at a finite 
h > 0. This will make us able to show that the return to the simplest VT — symmetric 
discrete models finds its unexpected reward in a complete suppression and elimination 
of the "unknown" annihilation and creation processes. In the other words we shall 
reinstall a firmer parallel between a simplifying role of VT— symmetry in both the 
bound-state and scattering-state hypothetical experimental arrangements. 

2 Solvable discrete models of scattering 

Let us consider the Hamiltonian H = H^ M \g) = Hq + W(g) of the doubly infinite 
matrix form where the non-vanishing part of the matrix W(g) = gV^ will be 



linear in the real coupling g and where the matrix V^ M > itself will be tridiagonal and 
formed just by the four off-diagonal nonvanishing matrix elements. These elements 
will be arranged in such a way that using the definition ([]]), the VT— symmetry of 
the complete Hamiltonian will be guaranteed, 

y(M) _ y(M) _ 1 y(M) _ y(M) ^ 

V l-M-M — V M-1,M ~ 1 ' V -M,l-M — v M,M-± ~ \ ' 

The resulting Hamiltonian H can be interpreted as a discrete kinetic-energy operator 
complemented by an interaction mimicking the VT— symmetrized pair of delta func- 
tions [9J. At the smallest "distances" M = 1, 2, . . . our model (jlj) may also resemble 
certain solvable short-range square- well differential-operator Hamiltonians [10]. In 
the free- motion case the above-mentioned connection between our H(0) = Ho and 
the Runge-Kutta Laplacean may be recalled to explain the origin of the constraint 
E G (0,4/h 2 ). This is a peculiarity which is well known in the bound-state context 
[2] . Here this restriction proves equally important for the physical consistency of the 
scattering boundary conditions (jSJ). 

In what follows, we intend to search for the solutions of Schrodinger eq. (121 + (!3]) 
using the standard matching method. We should emphasize that in the scattering 
scenario the key specific feature of wave functions is that they are constructed at 
any energy (from the allowed interval with, say, (p G (0, 7r)) and that they are not 
VT— symmetric themselves (this symmetry is broken by the boundary conditions). 
At the same time, due to the compact nature of the range of our interactions W, 
the non-compact character of the wave functions is fully characterized by eq. ([3]). 



Thus, in place of the doubly infinite matrix H(x) with the structure 





-1 








-1 


— 1 — X 

-l+x 


-1 








-1 




-1 








-1 


-l+x 

— 1 — X 


-1 



we only have to study the "central" submatrices of H in which W ^ 0. 

In principle, we could consider both the even- and odd- dimensional Ws. Never- 
theless, in the context of bound states we already saw that the difference between 
the 2M— and 2M + 1— dimensional cases is purely formal [UJ. For this reason we 
shall work just with odd dimensions here. This choice has the two marginal formal 
merits in containing the "first nontrivial" three-dimensional model at M = 1, 

yd) = 



o-io 

1 1 
0-10 



and in allowing the perceivably less puzzling indexing of the matrix elements by the 
parity-symmetric integers k = . . . , —2, —1, 0, 1,2,.. .. 



2.1 M = l 

At M — 1 the set of matching conditions involves just the following three rows of 
the central subset of the complete Schrodinger equation Htp = Eip, 



-1 2cos(/9 — 1 — x 
-1 + x 2 cos ip -1+x 
— 1 — x 2 cosy? —1 



Tew 
Te 2 ^ 



0. 



From the first and third row we get 1 + R = (1 + x)ip = T so that the remaining 
row multiplied by 1 + x, viz, equation 

(x 2 - 1) (e - ^ - + T e^) + 2T cos y? + (a; 2 - 1) T = . 



leads to the solution in closed form, 



1 



1 + L4' 



R 



-LA 



1 + iA' 



A 



1 - x< 



cot 



We may immediately verify that 



\R\ 2 + \T\ 2 = 1. 



This enables us to conclude that in spite of its non-Hermiticity, our scattering model 
conserves the probability at M — 1. 



2.2 M = 2 



At the next integer index M = 2 the set of matching conditions comprises the 
following five items, 



-1 2cosy? -l-x 

-l+x 2 cosy? -10 

-1 2cosv? -10 

-1 2cosy? -l + x 

-l-x 2cosy? -1 



e'^ + Re 3 ^ 
e~'^ + Re 1 * + x-i 

Te 3i V 



From the first and last line we get 



(l + x) X -i = -x{e- iv + Re i(p ), (l+x) X i = ~xTe ilp . 
This enables us to consider just the three modified matching conditions 

-l+x 2 2cosy? -10 
-1 2 cosy? -1 

-1 2cosy? -l+x 2 



e'^ + Re^ 
(l + x) ifj 
Tew 
Te 2 ^ 



0. 



The first row gives 



(1 + x )^ = l+ x 2 e- 2iv + (1 + x 2 e 2iLp ) R 



while the third row offers 



[l + x)4> = (l + x 2 e 2iLp )T 



so that we may eliminate ip and obtain the first rule for R and T, 

x 2 sin 2tp 



T = R+- —— = R+- _, A 



l + x 2 e 2 ^ 



1 + iA ' 



l+x 2 cos 2tp 



The remaining middle row leads to the third independent formula for 

1 + (R + T)e^ 

We may combine all three representations of (1 + x)ip and extract the second rule 
for R and T. In the light of the above representation of the difference T — i? we shall 
complement it by the second rule which determines the sum R + T. Such a recipe 
leads to the particularly compact final result, 

1 — ia 1 — \(5 



2R = 



1 + ia 1 + i/3 ' 



where 



m 1 - ia 1-iB 

2 T = I — 

1 + ia 1 + i/3 ' 



cos 2(/9 cot ip sin 2<^ 

a = 7^ P 



1 — 2x 2 cos 2 ip ' 1 + x 2 cos 2(/9 

Since both a; and /3 are real, it is immediate to prove that 

|i?| 2 + |T| 2 = 1. 

We see that in the model with M = 2 the flow of probability is conserved as well. 
One feels tempted to expect such a unitary-type behavior of the amplitudes at all 
the integer "interaction distances" M. 

Let us test such a conjecture on the next version of our model. 



2.3 M = 3 



Let us abbreviate U- m = e" miv + iZe™^ and L n = T e nilp and partition the seven 
matching conditions at M = 3 as follows, 









U-3 




C/_4 


2 cos y? — 1 — x 








U-2 + X-2 







— 1 + X 2 COSy? 


-1 






U-i + x-i 















^0 









-1 


2 cosy? — 1 + x 




£i +Xi 











— 1 — x 2 cos <f 




L 2 + X2 













U 




u 



The first and last lines give 

(1 + x)x-2 = -a;f7_2 , (1 + x)x 2 = -2^2 

and the elimination of the left-hand-side expressions gives the following reduced set 
of the five matching conditions, 



2 cosy? —1 




U-2 




(1 " X 2 )C/_ 3 


— 1 2 cos y? —1 




(1 +*)([/_! + X -l) 







— 1 2cosy? —1 




(l + x)^o 







— 1 2cosy? —1 




(l + x)(L 1 + xi) 







— 1 2 cosy? 




L 2 




(l-x 2 )L 3 



From the first and last equation we eliminate 

(1 + x) x-i = —xU-i + x 2 U-s , (1 + x)xi = —xLi + x 2 L 3 



and insert these expressions in the remaining three equations, with the result 



-1 2cosy? -10 
-1 2 cosy? -1 
0-1 2 cosy? -1 



U-2 
[/_! + X 2 U_ 3 
(l+x)l/>Q 

L x + x 2 L 3 
L 2 



Let us rewrite these equations again as the three non-equivalent definitions of 



(1 +x)4> = U + 2x 2 cosy?£/_ 3 , 



(1 + x)ifj = L + 2 x 2 cos y? L 3 , 



'l + x)ip = 



L x + x 2 L 3 + C/_! + C/_ 3 



2 COSy? 

and eliminate -?/>o m two alternative ways which define the difference 

1 + 2 x 2 e" 3 ^ cos if 1 - i7 



T-R = 



l + 2x 2 e 3i ^ cos 1 + i7 ' 



and the sum 



T + R = -e 



-2i v 1 - eiip cos V 5 - ^e -2 ^ cos 2ip 



1 — e~ il ^ cos ip — x 2 e 2lifi cos 2y? 
From these formulae it is again easy to derive 



l-Rl + \T\ 



i.e., the desirable conservation-of-probability law at M — 3. 



2.4 M = 4 

Out of the nine lines of the M = 4 matching conditions 
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2 cos y? — 1 — x 








+ X-2 







— 1 + X 2 COSy? 


-1 






+ X-i 















^0 









-1 


2 cosy? — 1 + x 
— 1 — x 2 cos y? 




L 2 + X2 














Lz + X3 

















L 5 



we may eliminate the first and last line using the general formula 

(1 + x)xi-m = -xUx-m , (1 + x)xm-i = -xL M -i . 

Also the rest of the solution can be perceived as a guide to the construction of the 
amplitudes R and T at any higher M. Indeed, once we return to the remaining seven 
matching conditions at M — 4, 



2 cosy? —1 
— 1 2 cosy? —1 



— 1 2 cosy? —1 
—1 2 cosy? 



U-3 

[l+x)(U- 2 + X-2) 

;i + *)([/_! + x_0 

(1 + x)ij 

(l + x)(L 1 + X i) 
(l + x)(L 2 + X 2) 
L 3 



(1- 


- x 2 )U_ A 



























(1 


-x 2 )L 4 



we may repeat the algorithm and eliminate its first and last line. Another general 
pair of formulae serves the purpose, 

(1 + x) X2-M = -x U 2 ~ M + x 2 U- M , (1 + x) XM~2 = -x L M _ 2 + x 2 L M 

after one inserts M — 4. In the subsequent step of the reduction procedure we arrive 
at the quintuplet of equations 



2 cosy? —1 




U-2 + X-2 




U- 3 /(l + x) 


— 1 2 cosy? —1 




U-i + x-i 







— 1 2cosy? —1 











— 1 2cosy? —1 




Li + Xi 







— 1 2 cosy? 




L 2 + X2 




^ L 3 /(l + x) 



Using the first and fifth equation again, we specify the last auxiliary quantities. 



( 1 + x) X-i = —x U-i + 2x 2 cos ip U-4 , ( 1 + x) xi — — x L\ + 2x 2 cos ip L A . 

This exemplifies the last step of the generic recurrent recipe because the next step 
will already involve the exceptional central element ipo- Thus, our knowledge of the 
expressions for x±i leads to the final triplet of conditions 

(1 + x)ip = U + x 2 (l + 2 cos 2ip)U- 4 , 
(1 + x)ip = L + x 2 (l + 2 cos 2f)L A , 

(1 + x)^ = Ll + U ~ 1 + x 2 (L 4 + U-t) . 
2 cosy? 

After the two alternative eliminations of ipo we routinely arrive at our last two linear 
equations for the two unknown quantities R+T and T — R. Their elementary though 
a bit clumsy solution will not be displayed here anymore. Whenever asked for, the 
proof of the conservation law at M = 4 as well as the further, more or less routine 
though increasingly tedious continuation of our construction to the higher "distances 
M between interactions" are left to the readers. 



3 Summary 



One of the most pleasant and encouraging observations made during many practical 
applications of quantum theory is that our basic understanding of experimental data 
can often be provided by fairly elementary mathematical models. Among them, a 
prominent role is played by the one-dimensional Schrodinger equation. Of course, 
the detailed physical interpretation of such a class of models can vary with the 
experimental setup and may range from the naive fitting scenario up to a schematic 
reduction of field theory to zero dimensions. 

In the latter, highly speculative context Bender and Milton [12] and Bender and 
Boettcher [13] revealed that phenomenological as well as theoretical purposes could 
be served very well by complex potentials exemplified by V(x) = ix 3 and supporting 
real spectra of bound states [H] . Later on, it has been clarified that the transition to 
the complex V(x) does not in fact violate any rules of Quantum Mechanics because 
even for complex potentials the Hamiltonian can be reinterpreted as self-adjoint after 
a suitable adaptation of the Hilbert space of states [15]. 

Jones [B] was probably the first author who analyzed the possibilities of the 
same adaptation of the Hilbert space in the scattering scenario. Although he choose 
one of the simplest and best understood potentials, viz., the delta function with a 
complex coupling, his conclusions concerning both the mathematical feasibility and 
the physical clarity of the complexified scattering problem were rather discouraging. 
His construction revealed that in spite of the ultralocal form of his toy model the 
scattered waves proved perceivably and counterintuitively distorted. 

In our present note we reanalyzed the situation by incorporating, in explicit 
manner, the postulate of the so called VT— symmetry of the Hamiltonian which is 
often being implemented in the constructive description of bound states in unusual 
Hilbert spaces. For this purpose we introduced and solved and entirely new class 



of discrete models of scattering. We were really surprised when we revealed that 
these models behaved differently in comparison with their similar VT— asymmetric 
predecessors of refs. [HI [7]. 

The key merit of our present family of models should be seen in the fact that 
not quite expectedly, they fully conserve the probability and do not seem to exhibit 
any signs of an asymptotic non-locality Moreover, since they are simple and exactly 
solvable, the emerging possibilities of their entirely standard practical applications 
and/or theoretical probabilistic interpretation do not seem to be an artifact of their 
present discretized mathematical form. 

We believe that on the background of certain pessimistic physics-related perspec- 
tives as formulated in refs. [51 E], our present results could serve as a source of new 
optimism, needed for the continuation of the search for some new manifestly non- 
Hermitian models of scattering. One can hope that the user-friendly features of our 
models will survive their extensions, both in the sense of returning to the continu- 
ous limit h — > and in the sense of finding their more-parametric descendants of a 
greater descriptive flexibility. 
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